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According to recent general-relativistic simulations, the coalescence of two spinning black holes 
(BHs) could lead to recoil speeds of the BH remnant of up to thousands of km s^" 1 as a result of 
the anisotropic emission of gravitational radiation. Such speeds would enable the merger product to 
escape its host galaxy. Here we examine the circumstances resulting from a gas-rich galaxy merger 
under which the ejected BH would carry an accretion disk with it and be observable. As the initial 
BH binary emits gravitational radiation and its orbit tightens, a hole is opened around it in the 
disk which delays the consumption of gas prior to the eventual BH ejection. The punctured disk 
remains bound to the ejected BH within the region where the gas orbital velocity is larger than the 
ejection speed. For a ~ lO 7 A/0 BH the ejected disk has a characteristic size of tens of thousands 
of Schwarzschild radii and an accretion lifetime of ~ 10 7 years. During that time, the ejected BH 
could traverse a considerable distance and appear as an off-center quasar with a feedback trail along 
the path it left behind. A small fraction of all quasars could be associated with an escaping BH. 



Introduction. The past few years witnessed a break- 
through in numerical relativity as simulations were able 
to follow for the first time the final coalescence phase of 
a binary black hole (BH) system due to the emission of 



gravitational radiation [1|, |2j, H, 0, H, 0] • While the re- 
coil speed of the remnant from binaries of non-spinning 
BHs is modest (< 200 km s -1 ), the coalescence of spin- 
ning BHs of nearly equal masses could yield recoil speeds 
of up to thousands of km s _1 0,0,0- The largest speeds 
are obtained for special orientations of the spins relative 
to the orbital plane, and it is unclear whether they occur 
in nature Q. To find out if they do, it is necessary to 
identify observational signatures of the ejected BHs. 

Accretion of interstellar medium (ISM) gas by a fast- 
moving BH produces only faint luminosities. For a 
BH ejection speed i> C j which is well above the sound 
speed of the ISM, the Bondi [9( accretion rate is Mb = 
7 x 10^ 7 n M^Vg 3 M Q yr" 1 , where M 7 is the BH mass 
in units of 10 7 M Q , vg is the BH speed in units of 
10 8 cm s _1 = 10 3 kms -1 , and n is the ISM den- 
sity in units of 1 cm~ 3 . Even if a thin accretion disk 
forms around the ejected BH and the Bondi accretion 
rate is converted to radiation with a high efficiency 
e = 0.1e_i, the resulting luminosity L = eMbc 2 = 
4x 10 39 e_in M|u^ 3 erg s _1 , will be difficult to detect at 
cosmological distances, except on the rare occasion when 
the ejected BH passes through a dense molecular cloud. 

The ejected BH could appear much brighter if it car- 
ries an accretion disk with it. As long as the disk mass 
is much smaller than the BH mass and the gas within 
the disk is orbiting at a speed far greater than the BH 
ejection speed (which is possible since v e j <C c), the gas 
will preserve the adiabatic invariants of its orbit around 
the BH and follow the BH along its ejection trajectory. 

Next we characterize the properties of accretion disks 
that could be carried by ejected BHs, and in the final 
section we discuss their observational signatures. 



Disk Parameters. Hydrodynamic simulations indi- 
cate that a gas-rich merger between a pair of galaxies 
tends to drive their supermassive BHs together with some 
gas towards the center of the merger product 13, II , 10] • 
The latest simulations indicate that the cold gas is not 
expelled by the feedback from early accretion episodes 
when the two BHs are far apart, and so a gaseous enve- 
lope forms around the coalescing BH binary. 

The coalescence of the two BHs is driven at first by dy- 
namical friction on the gas and stars 11, 12, 13| • Once 
the BH binary separation shrinks to the regime where its 
orbital speed exceeds ~ 10~ 2 c, gravitational radiation 
alone could cause coalescence within a Hubble time 14 1 . 



Although the binary might stall at larger separations due 
to the emptying of the loss cone of stars that could ex- 
tract angular momentum from its orbit (the so-called fi- 
nal par sec problem), there are many plausible processes 
that would refill this loss cone [3] or enable gas to grind 
the orbit down to the gravitational-radiation dominated 
The end result is likely to be a binary BH 
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regime 

system that continues to tighten on its own due to the 
emission of gravitational radiation, surrounded by an en- 
velope of gas [iH . The gas will likely settle to a cold disk 
around the BH binary, since its cooling time is short. The 
disk may include some gas that was originally associated 
with either BH individually but joined into a common 
envelope once the binary separation was sufficiently re- 
duced. 

Recent hydrodynamic simulations of a common co- 
rotating disk of gas around a binary system of BHs of 
comparable mass [HI [l6| indicate that a hole tends to 
open across a region of a radius equal to twice the bi- 
nary semi-major axis a. The clearing of the hole is sim- 
ilar to the opening of a ga p b y a massive planet in a 
gaseous disk around a star [ljj. The torque exerted by 
the binary BH on the outer edge of the hole, pushes gas 
elements outwards at that location 
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Hil ]. The cor- 
responding transfer of angular momentum promotes the 
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binary coalescence process 16, [l^|. As the binary sep- 
aration shrinks, the hole is expected to gradually close 
from the outside due to the viscous transport of angular 
momentum by the _gas. The viscous timescale of the disk 
can be expressed 21 1 in terms of the a-parameter , 



£ visc (r) =4.1 x 10 4 
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Mi 



yr, 
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where M 7 is the total BH binary mass in units of 10 7 M Q , 
r3 is the radius r relative to the binary center of mass 
in units of 10 3 Schwarzschild radii of M (= 2GM/c 2 ), 
a_i is the viscosity parameter scaled to a characteris- 
tic value of 0.1, the opacity is assumed to be dominated 
by Thomson scattering for a primordial gas composition, 
and r\ = (e/0.1) /(L/Le) with e being the radiative ef- 
ficiency and L being the disk luminosity in units of the 
Eddington limit L E = 1.4 x 10 45 M 7 erg s" 1 if the a-disk 
were to extend down to the last stable orbit around a 
single BH. For comparison, the time it takes a binary on 
a circular orbit to coalesce due to the emission of gravi- 
tational radiation is [201 



2xl0 6 f|Qa 4 M 7 yr, 
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where [i = (M1M2/M) is the reduced mass, and M = 
(Mi + M2) is the total mass of the binary with BH 
masses Mi and M2. As before, a subscript (...)3 denotes 
a lengthscale in units of 10 3 Schwarzschild radii of M. 
Hereafter, we focus our attention on the case where the 
BH masses are comparable (M/fi ~ 4) since the recoil 
speed is small when one of the binary members is much 
lighter than the other pj. 

There is a minimum radius to the cavity around the 
binary, i?j n , below which the decay time of the binary 
orbit will be much shorter than the viscous time required 
to refill the cavity [ill (23[. We obtain this radius by 
requiring i gw < i v i sc and substituting r ~ 2a for the 
radius of the hole around the BH binary [3]. This gives, 



Rin,3 ~ 0.65M£ 



0.077, 
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As the binary orbit continues to tighten with a < 
0.5i?i n , the disk inner radius does not have sufficient time 
to close much farther before the binary BH coalesces. 
The ejected BH would therefore carry an initially punc- 
tured disk with an inner cavity radius of i?;,, and an outer 
radius of i? Q ut 



GM/Ugj or equivalently [34 1 
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We therefore find that a bound disk could survive (i.e. 
i?ou t > i?in) for the relevant regime of recoil speeds 
[35| . Well interior to R ou t, the orbital velocity of the 
gas exceeds v C j and the gas elements maintain the adia- 
batic invariants of their orbit around the relatively "slow- 
moving" BH. This would hold even if the BH is kicked in 
the plane of the disk. The gas outside i? ou t is unable to 



respond sufficiently quickly to the ejection of the BH and 
so it is left behind. After the BH remnant is ejected, the 
captured disk relaxes to a viscous equilibrium state and 
fills its central cavity on a timescale i V i S c(^m), which is 
much shorter than its global accretion lifetime t v i sc (R out ). 
Magnetohydrodynamic (MHD) simulations of an initially 
toroidal gas distribution demonstrate this behaviour (see, 
e.g. Ref. [25| and references therein). 



The total mass of the a-disk [2l| around the ejected 
BH is, 



M disk « 1.9 x 10 6 a:;-V°- 6 M 7 2 - 2 ^ 
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The above expression is valid only in the regime where 
-Wdisk "C M since we ignored the mass of the disk (as well 
as its self-gravity) in the overall momentum balance of 
the BH+disk system. The condition that the BH would 
not carry a disk mass in excess of M implies 

i?out,3 < flmax,3 = 1.5 X lOVf V ~ 0A3 M^ 86 . (6) 

For the characteristic surface density and pressure of an 
a-disk [2lJ, we find that the ram-pressure of the ISM 
through which it passes (as well as the amount of ISM 
mass intercepted) can be ignored. 

Observable Signatures. A gas-rich merger of two 
galaxies naturally creates an environment that includes 
a binary BH system surrounded by gas near the center 
of the merger product 11 , IjJ ■ As gas accumulates out- 
side the binary, it is likely to form a punctured disk that 
would survive through the binary coalescence process. 

The lifetime of a disk with an outer radius i? ut given 
by Eq. (4} can be derived ^ from Eq. flU, 



t disk « 8.4 x 10 6 aZ°i 8 v° A M}- 2 v 8 
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yr. 



(7) 



During the disk lifetime tdisk, the ejected BH could tra- 
verse a distance of 



vejidisk « 8.6 aI?-y- 4 M 



(8) 



and appear as an off-center quasar. An offset of d ~ 10 
kpc can be resolved, as it corresponds to angular scales of 
order an arcsecond at cosmological distances. A notice- 
able displacement might exist even for BHs which remain 
bound to their host galaxy halo. There has already been 
a claim for a detection of a displaced quasar [27|], which 
has been disputed by subsequent analysis [11] . An undis- 
puted example for the preceding phase of a compact BH 
binary was reported recently in a different system [29| . 

Interestingly, the offset in Eq. ([5]) increases as the 
ejection speed decreases. This scaling originates from 
the steep inverse dependence of the disk size on the 
ejection speed. However, it is only valid as long as 
v C j ^> 200 km s _1 , because the accretion disk is not ex- 
pected to follow the ejected BH from radii where the gas 
orbital speed is lower than ~ 200 km s _1 at which the 
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2oo is typically 
and hence is 



tidal force from the galaxy is important. For low ejec- 
tion speeds, v e j < a, the BH would only carry the frac- 
tion of the disk interior to i? out ~ mm{(GM/a 2 ), -R max }, 
where a is the velocity dispersion of the host galactic 
spheroid. Expressing a = 200(T2oo km s -1 , the value of 
(GM/ct 2 ) = 1.1 x 10 3 x (10 3 x 2GM/c 2 )a. 
larger than i? max for M 7 > 1 in Eq. 
of secondary importance. At low ejection speeds, the 
deceleration of the BH by the external gravitational po- 
tential of the galaxy and by dynamical friction on the 
background stars and gas, would limit the displacement 
of the slow-moving BH and tend to bring it back to the 
center of its host galaxy within a few dynamical times 

For a high ejection speed, v$ > 1, the displacement of 
the ejected BH relative to the host galaxy would also be 
noticeable in velocity (redshift) space. The velocity offset 
could be detected if the material bound to the ejected 
BH (including gas and possibly some captured stars) is 
capable of producing emission lines. 

The typical lifetime of quasars is estimated t 31] to be 
comparable to the value of tdisk in Eq. ijTj). However, for 
equal-mass binaries with random spins (uniformly dis- 
tributed between dimensionless values of and 0.9 and 
of random orientations) only ~ 15 ± 6% of all mergers 
are expected to result in an ejection speed vs > 0.5 and 
only - 3 ± 2% in v s > 1 Q, above the escape speed of 
a massive galaxy. Hence, we expect only a small frac- 
tion of all quasars to be associated with an unbound 
BH that escapes the galaxy. Indeed, the fact that al- 
most all nearby galaxies possess a nuclear BH [32I even 
though some formed out of gas-poor mergers (from which 
a new BH could not have been created to replace an ex- 
pelled BH) implies that unbinding kicks are rare. Simi- 
larly to most bright quasars, the luminosity of a kicked 
BH could be close to the Eddington limit, except that 
its fuel reservoir is short-lived and cannot be replenished 
unless v C j < a and the BH settles back to the center of 
the galaxy. 

The motion of the BH and its displacement to regions 
of low gas density would reduce the growth rate of its 
mass relative to the case where the remnant BH stays 
at the center of its host galaxy. By the time an ejected 
BH is driven back to the center by dynamical friction, a 
substantial fraction of the cold gas there may already be 
converted into stars or expelled in a galactic wind. 

The feedback on the host galaxy from the radiation or 
wind produced by the accreting disk around the ejected 
BH would obviously depend on the ejection speed. Pre- 



vious treatments of BH feedback in galaxy mergers [11 1 
ignored the gravitational-wave recoil and assumed that 
the BH remnant stays fixed at the center of mass of 
its host galaxy. The energy/momentum feedback was 
therefore confined to the innermost (and hence densest) 
resolution element of gas. However, the actual proba- 
bility distribution of ejection speeds Q allows for sub- 



stantial BH displacements in many cases, for which the 
energy /momentum output from the BH would interact 
differently with the surrounding gas. If the ejected BH 
carries an accretion disk as discussed in this paper, its 
feedback on the host galaxy would be different than pre- 
viously calculated. The distributed energy/momentum 
release of a displaced BH would leave a trail of evidence 
that traces its origin to the center of the galaxy through 
the mark of its feedback along its path. For example, 
this mark could be traced by observing enhanced Ha line 
emission due to the ionization trail imprinted by the BH 
on the diffuse ISM along its trajectory. 

The basic conclusions of this Letter would be altered 
if gas is expelled from the vicinity of the BH binary by a 
powerful galactic wind, driven by supernovae or quasar 
activity prior to the binary merger event. A burst of star 
formation activity is generically expected to accompany 
the merger process [ll| and could affect the assembly of 
cold gas into the circumbinary accretion disk that was 
postulated in this Letter. It is also possible that low- 
level accretion across the gap surrounding the BH binary, 
would consume the gas in the initial disk and weaken the 
late accretion episode considered here. 

In the future, it would be useful to improve upon the 
approximate treatment presented here by simulating the 
dynamical response of an MHD disk to the final stage of 
the BH binary coalescence process and the subsequent 
recoil of the BH remnant. Such a calculation can be 
done with existing codes [llj [l6| that were used so far 
to address other aspects of the binary coalescence phe- 
nomenon. Cosmological hydrodynamics codes [H| could 
also be used to simulate the feedback trail that an ejected 
BH imprints on the ISM of its host galaxy, which is remi- 
niscent of the track of an elementary particle in a bubble 
chamber. 
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